Identification of Genome-Wide Variants and Discovery of Variants Associated with Brassica rapa Clubroot Resistance Gene Rcr1 through Bulked Segregant RNA Sequencing.

Clubroot, caused by Plasmodiophora brassicae, is an important disease on Brassica species worldwide. A clubroot resistance gene, Rcr1, with efficacy against pathotype 3 of P. brassicae, was previously mapped to chromosome A03 of B. rapa in pak choy cultivar "Flower Nabana". In the current study, resistance to pathotypes 2, 5 and 6 was shown to be associated with Rcr1 region on chromosome A03. Bulked segregant RNA sequencing was performed and short read sequences were assembled into 10 chromosomes of the B. rapa reference genome v1.5. For the resistant (R) bulks, a total of 351.8 million (M) sequences, 30,836.5 million bases (Mb) in length, produced 120-fold coverage of the reference genome. For the susceptible (S) bulks, 322.9 M sequences, 28,216.6 Mb in length, produced 109-fold coverage. In total, 776.2 K single nucleotide polymorphisms (SNPs) and 122.2 K insertion / deletion (InDels) in R bulks and 762.8 K SNPs and 118.7 K InDels in S bulks were identified; each chromosome had about 87% SNPs and 13% InDels, with 78% monomorphic and 22% polymorphic variants between the R and S bulks. Polymorphic variants on each chromosome were usually below 23%, but made up 34% of the variants on chromosome A03. There were 35 genes annotated in the Rcr1 target region and variants were identified in 21 genes. The numbers of poly variants differed significantly among the genes. Four out of them encode Toll-Interleukin-1 receptor / nucleotide-binding site / leucine-rich-repeat proteins; Bra019409 and Bra019410 harbored the higher numbers of polymorphic variants, which indicates that they are more likely candidates of Rcr1. Fourteen SNP markers in the target region were genotyped using the Kompetitive Allele Specific PCR method and were confirmed to associate with Rcr1. Selected SNP markers were analyzed with 26 recombinants obtained from a segregating population consisting of 1587 plants, indicating that they were completely linked to Rcr1. Nine SNP markers were used for marker-assisted introgression of Rcr1 into B. napus canola from B. rapa, with 100% accuracy in this study.

Fine mapping of Rcr1 and analyses of its effect on transcriptome patterns during infection by Plasmodiophora brassicae.

BACKGROUND: The protist Plasmodiophora brassicae is a biotrophic soil-borne pathogen that causes clubroot on Brassica crops worldwide. Clubroot disease is a serious threat to the 8 M ha of canola (Brassica napus) grown annually in western Canada. While host resistance is the key to clubroot management, sources of resistance are limited. RESULTS: To identify new sources of clubroot resistance (CR), we fine mapped a CR gene (Rcr1) from B. rapa ssp. chinensis to the region between 24.26 Mb and 24.50 Mb on the linkage group A03, with several closely linked markers identified. Transcriptome analysis was conducted using RNA sequencing on a segregating F1 population inoculated with P. brassicae, with 2,212 differentially expressed genes (DEGs) identified between plants carrying and not carrying Rcr1. Functional annotation of these DEGs showed that several defense-related biological processes, including signaling and metabolism of jasmonate and ethylene, defensive deposition of callose and biosynthesis of indole-containing compounds, were up-regulated significantly in plants carrying Rcr1 while genes involved in salicylic acid metabolic and signaling pathways were generally not elevated. Several DEGs involved in metabolism potentially related to clubroot symptom development, including auxin biosynthesis and cell growth/development, showed significantly lower expression in plants carrying Rcr1. CONCLUSION: The CR gene Rcr1 and closely linked markers will be highly useful for breeding new resistant canola cultivars. The identification of DEGs between inoculated plants carrying and not carrying Rcr1 is an important step towards understanding of specific metabolic/signaling pathways in clubroot resistance mediated by Rcr1. This information may help judicious use of CR genes with complementary resistance mechanisms for durable clubroot resistance.

Heteroconium chaetospira induces resistance to clubroot via upregulation of host genes involved in jasmonic acid, ethylene, and auxin biosynthesis.

An endophytic fungus, Heteroconium chaetospira isolate BC2HB1 (Hc), suppressed clubroot (Plasmodiophora brassicae -Pb) on canola in growth-cabinet trials. Confocal microscopy demonstrated that Hc penetrated canola roots and colonized cortical tissues. Based on qPCR analysis, the amount of Hc DNA found in canola roots at 14 days after treatment was negatively correlated (r = 0.92, P<0.001) with the severity of clubroot at 5 weeks after treatment at a low (2×10(5) spores pot(-1)) but not high (2×10(5) spores pot(-1)) dose of pathogen inoculum. Transcript levels of nine B. napus (Bn) genes in roots treated with Hc plus Pb, Pb alone and a nontreated control were analyzed using qPCR supplemented with biochemical analysis for the activity of phenylalanine ammonia lyases (PAL). These genes encode enzymes involved in several biosynthetic pathways related potentially to plant defence. Hc plus Pb increased the activity of PAL but not that of the other two genes (BnCCR and BnOPCL) involved also in phenylpropanoid biosynthesis, relative to Pb inoculation alone. In contrast, expression of several genes involved in the jasmonic acid (BnOPR2), ethylene (BnACO), auxin (BnAAO1), and PR-2 protein (BnPR-2) biosynthesis were upregulated by 63, 48, 3, and 3 fold, respectively, by Hc plus Pb over Pb alone. This indicates that these genes may be involved in inducing resistance in canola by Hc against clubroot. The upregulation of BnAAO1 appears to be related to both pathogenesis of clubroot and induced defence mechanisms in canola roots. This is the first report on regulation of specific host genes involved in induced plant resistance by a non-mycorrhizal endophyte.

Patterns of differential gene expression in Brassica napus cultivars infected with Sclerotinia sclerotiorum.

SUMMARY The fungal pathogen Sclerotinia sclerotiorum infects a broad range of dicotyledonous plant species and causes stem rot in Brassica napus. To elucidate the mechanisms underlying the defence response, the patterns of gene expression in the partially resistant B. napus cultivar ZhongYou 821 (ZY821) and the susceptible cultivar Westar were studied using a B. napus oligonucleotide microarray. Although maximum differential gene expression was observed at 48 h post-inoculation (hpi) in both cultivars, increased transcript levels were detected in cv. ZY821 at the earlier stages of infection (6-12 hpi) for many genes, including those encoding defence-associated proteins, such as chitinases, glucanases, osmotins and lectins, as well as genes encoding transcription factors belonging to the zinc finger, WRKY, APETALA2 (AP2) and MYB classes. In both cultivars, genes encoding enzymes involved in jasmonic acid, ethylene and auxin synthesis were induced, as were those for gibberellin degradation. In addition, changes in the expression of genes encoding enzymes involved in carbohydrate and energy metabolism appeared to be directed towards shuttling carbon reserves to the tricarboxylic acid cycle and generating reactive oxygen species. Transcripts from genes encoding enzymes involved in glucosinolate and phenylpropanoid biosynthesis were highly elevated in both cultivars, suggesting that secondary metabolites are also components of the response to S. sclerotiorum in B. napus.